TT D —— — a it — A LE 


LL Á- -_—_—_- 
— S 


Michael J. Donoghue,’ 
Richard G. Olmstead,’ 
James F. Smith,’ and 
Jeffrey D. Palmer’ 


PHYLOGENETIC 
RELATIONSHIPS OF 
DIPSACALES BASED ON 
rbcL SEQUENCES! 


ÅBSTRACT 


T pele sequences of the chloroplast gene rbcL were obtained from Lonicera, Sambucus, Adoxa, and Cornus. 
Á ylogenetic analyses of these sequences, along with previously sequenced Dipsacales and placeholders for other 
eee and Rosidae, lead to the following conclusions: (1) the genera of Caprifoliaceae (in any traditional sense) 
© not form a monophyletic group; (2) Symphoricarpos and Lonicera (representing Caprifoliaceae sens. str.) are 
related and are united with Valerianaceae and Dipsacaceae; and (3) Adoxa and Sambucus are directly linked and 
are possibly related to Viburnum. Our analyses also cast doubt on the monophyly of Dipsacales and Asteridae. These 


findings are generally congruent with other molecular results and with hypotheses based on morphological data. 


Uncertainty has long surrounded the phyloge- 
netic relationships of taxa traditionally assigned to 
the angiosperm order Dipsacales (Donoghue, 1983). 
It has been unclear, for example, whether Adoxa 
is related to Dipsacales, whether Sambucus and 
Viburnum are closely related to one another or to 
genera of Caprifoliaceae sens. str., or whether some 
Caprifoliaceae are more closely related to Valeri- 
anaceae and/or Dipsacaceae than they are to each 
other. The broader phylogenetic relationships of 
these plants have also been debated. It is still un- 
certain whether Dipsacales were derived from (i.e., 
Ng nested within) Rubiales, as implied by Cronquist 
981, 1988), or whether they (and possibly As- 
erales) are more closely related to rosid groups 
such as Cornales, as suggested recently by Takh- 
tajan (1987; also see Dahlgren, 1983; Thorne, 
1983), 

‘a Donoghue (1983) considered morphological ev- 

ence bearing on these questions, especially the 
relationships of Viburnum. However, owing to un- 
atria over broader phylogenetic relationships, 
© was able to reach only tentative conclusions, 


conditional on particular assumptions. For exam- 
ple, if Caprifoliaceae sens. lat. and Dipsacales are 
both monophyletic, then morphological data di- 
rectly unite Viburnum with Sambucus plus Adoxa, 
whereas under some other assumptions (e.g., Dip- 
sacales polyphletic) this connection is unclear. Sub- 
stitution of different outgroup arrangements (sensu 
Donoghue & Cantino, 1984) yielded few unequiv- 
ocal conclusions. Happily, however, the position of 
Adoxa was an exception. If Adoxa is assumed to 
be related in some way to other Dipsacales (as 
opposed, for example, to Saxifragaceae), then it is 
invariably linked with Sambucus on the basis of 
several morphological features, including Adoxa- 
type embryo sac development. 

Here we present a preliminary analysis of the 
phylogenetic relationships of dipsacalean taxa based 
on molecular data. In particular, we report the 
results of parsimony analyses of chloroplast rbcL 
sequences (see Ritland & Clegg, 1987; Palmer et 
al., 1988; Soltis et al., 1990; Doebley et al., 1990; 
Olmstead et al., 1992) of various Dipsacales and 
possibly related groups, including new sequences 
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from representatives of four critical genera: Adoxa, 
Sambucus, Lonicera, and Cornus. Although we 
are able to make some comparisons between the 
molecular results and those based on morphological 
characters, a detailed comparison must await the 
completion of an expanded morphological study 
(Donoghue, in prep.) and sequencing of additional 
Dipsacales (in progress). In the meantime, we high- 
light conclusions that are most strongly supported 
by rbcL data, so that these can be taken into 
account in formulating a phylogenetic system of 
Asteridae. 


MATERIALS AND METHODS 


rbcL sequences were newly obtained from four 
taxa: Adoxa moschatellina L. of Adoxaceae, Cor- 
nus mas L. of Cornaceae, and Sambucus racemosa 
subsp. sieboldiana (Miquel) Hara, and Lonicera 
prolifera Rehder of Caprifoliaceae sens. lat. The 
Wisconsin Department of Natural Resources kind- 
ly permitted the collection of Adoxa specimens in 
Vernon County, Wisconsin (Donoghue, Smith, & 
Kuetchenreuther 1, WISC). Material of the other 
species was obtained by M. Donoghue at the Arnold 
Arboretum, Jamaica Plain, Massachusetts (where 
vouchers are deposited). In each case, fresh leaves 
were collected on ice and then stored at —80°C, 
and DNA extractions were carried out using a 
modified CTAB method (Doyle & Doyle, 1987; 
Smith et al., 1992). 

A segment of the chloroplast genome containing 
the rbcL gene was amplified by the polymerase 
chain reaction (PCR), using the primers described 
in Olmstead et al. (1992; see their fig. 1). This 
procedure yielded double-stranded segments of ap- 
proximately 1550 bp. Each amplification product 
was then digested with the restriction endonucleas- 
es EcoRV and Xbal, known to recognize other 
Asteridae rbcL at nucleotide positions 103-108 
and 1269-1274, respectively. As expected, this 
procedure yielded three fragments in Adoxa, Lo- 
nicera, and Cornus: a 5' segment of 102 bp, an 
internal segment of 1172 bp including most of the 
rbcL gene, and a 3’ segment of approximately 290 
bp (minor variation owing to length difference 
downstream of rbcL). However, in the case of 
Sambucus the double digest yielded four fragments 
due to an additional EcoRV site present at nucle- 
otide positions 1018-1023 bp. This site is also 
found in Vibrunum plicatum, and more extensive 
restriction site studies indicate its presence in other 
(but not all) species of Viburnum (Donoghue & 
Sytsma, unpublished). 

The fragments obtained were cloned in the plas- 


mid vector BlueScript Sk+ (Stratagene, Inc.) as 
described by Olmstead et al. (1992), and single 
stranded DNAs were prepared (phenol-chloroform 
extraction) from selected colonies screened for in- 
serts of the larger EcoRV-Xbal fragments (or the 
equivalent fragments in Sambucus). Sequencing 
was carried out using the dideoxy chain termination 
technique, with internal primers kindly provided 
by G. Zurawski (DNAX Research Institute, Palo 
Alto, California). In some cases the smaller 5’ and 
3' fragments were recovered and sequenced, so 
that the entire rbcL sequence was obtained. How- 
ever, the analyses reported below are based solely 
on the internal 1172 bp segment. Using the known 
tobacco sequence as a reference, differences in our 
sequences were easily scored; length differences 
were never encountered in the region of interest. 
The sequences utilized in the present study will be 
available in GenBank, and may also be obtained 
directly from M. Donoghue or R. Olmstead. 

A total of twelve taxa were included in the phy- 
logenetic analyses reported here. These were se- 
lected so as to represent Dipsacales and other major 
lines of Asteridae and relevant Rosidae, as well as 
to include suggested relatives of the dipsacalean 
taxa, especially Adoxa. In addition to the four 
newly sequenced taxa discussed above, seven 5°- 
quences were obtained from workers in J. paion? 
laboratory: R. Olmstead provided Hedera helix L. 
(Araliaceae) and Symphoricarpos albus yi 
(Caprifoliaceae sens. str.) (see Olmstead et p 
1992), as well as Spigelia marilandica L. ( 
ganiaceae) and Menyanthes trifoliata L. eni 
anthaceae) (R. Olmstead, unpublished); and H. i 
Michaels (pers. comm.; Michaels et al., in A 
is responsible for Viburnum acerifolium L. (Cap 
rifoliaceae sens. lat.), Dipsacus sativus Honckeny 
(= D. fullonum L.) (Dipsacaceae), and Ve EA 
officinalis L. (Valerianaceae). The published A 
quence (Soltis et al., 1990) of Heuchera micran 
Dougl. (Saxifragaceae) was also included. es 

Elimination of invariant nucleotide positions 
194 variable sites, and by omission of passes 
phies this was further reduced to 88 poten y 
informative characters (Table 1). Of the latter, is 
are binary, 12 are unordered 3-state, and 2 . 
unordered 4-state characters, so there are 4 er 
imum of 104 character state changes (steps) in © 
absence of homoplasy. Note that the vast majority 
of potentially informative sites occur at piei 
sitions within codons; only eight are at first post 
and six at second positions. sk 

Character state changes were equally weigh i 
in most of our analyses, but we did experimen" 
with several differential weighting schemes, mc 
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Taste 1. Twelve taxa scored for 88 potentially informative rbcL characters (invariant sites and autapomorphies 
removed), with characters numbered by base position within the gene. All positions are third positions within codons, 
unless otherwise indicated above the character number. The complete set of potentially informative nucleotides is 
given for Heuchera (used to root trees); a nucleotide is shown for other taxa only if it differs from Heuchera (a dot 
indicates that the same nucleotide is present). Missing information (due to ambiguity in scoring) is indicated by ‘‘?.” 


CS ee eee 


nN Coena = SRF 
AAAA ITEE A AE A A 
e a a TT TTT E E E R OE a a 
HEUCHERA TAGGTTCCATCGAGTCCCGCCCTITGTCTTAACAACTTACCACAGAGA 
HEDERA ke. E . GC A A Bios: Č G IR ae A AN: 
CORNUS Tithe oe KES A Tet es G SMA re T T O e 
SPIGELIA C?.AC..TGCTCGAA oR. oO OG.) O eat hone eet ae 
MENYANTHES SAR REE: ssh hat Os ne eT R ce Coe ee ae G.G 
ADOXA ES GOT C2 Pe KR r T Rs, Dec: A N ren fh Oh Ge E A G 
SAMBUCUS ee ae ee eae ee Tree T ee War ta a 0 oe G 
VIBURNUM A NE AT AAO A M aE Oe Ve a DEA AN A I E 
LONICERA CG A ehe N ar Be Let De Gi. aaa G.TCAT. 
SYMPHORICARPOS , GA. . A BC NC GEROT- ENOC TOOG  OGA. 6. TEA. r @. , 
VALERIANA ike EA. OOS Thc. OT Te bee Oe eam Ga WS eek Seen RIM. 
DIPSACUS C. AACA OA = kA T teg AGGT..A G.TC rhe 
PESG Gece ee wee eRe leseesaser oes ene E 
KRRRRRRQDODDDDDDDBDOAMAAAADeRK KKK KrrKrKr rr rr rr Kerr ee 
HEUCHERA CACGGCCTTIGTGATTACCTTCGCATCTACGATTATATAAAT 
DEMA G O T OF Ge. Gs lane A EE AET N 
ONUS G. Ea: AaS ATAOE Ga. T Tae e GA 
SPIGELIA CANT CMI CA ONA aE Te. eT c.cG OSGi I 
MENYANTHES 6.2. .THeen..0@..c aah PA eh Ee ery, c.GC..G.CA 
ADOXA -GAA..CCT.¢cG.c hE tet te GAGC. GC _. GA 
SAMBUCUS CGAA..C. CG.cG..Gc.TTGCT..6.GC.6C CA 
VIBURNUM pet pio, oa OSES) a O O GOT. eRe OK. . CA 
LONICERA TCG. C..cCCaccG. .GC.A.G.T.&. S73 Cpe rọ 
mae T OG c ccGcce ERE E E a cet on L T AE 
ah T E AE Sete Gh GCTCGTAGCCG. CCGGTC 
T Ve sri et ryt Gia'¢ AKAN oe Gor. Go SUE Gs CC. eTO 


ing omission or down-weighting of various positions retention index (RI; Farris, 1989; also see Archie, 
within codons and weighting of transversions (Tv) 1989a). 
versus transitions (Ts). The Tv: Ts weighting of The lengths and implications of alternative tree 
1.3:1.0 used in the present study is based on topologies were explored using the Topological 
calculations by V. A. Albert (pers. comm.; see Constraints option in PAUP, and using finalized 
Albert & Mishler, in press), which in turn rely on portions of MacClade test version 2.99.11 (to be 
Substitution rates and Ty: Ts ratios for rbcL ob- published as MacClade 3.0; Maddison & Maddison, 
tained from a 5-taxon tree of land plants (Albert 1992). MacClade was also used to calculate char- 
et al., 1992). acter transformations of various types on trees of 
In view of the small number of taxa under con- interest (e.g., transversions versus transitions) and 
sideration, parsimony analyses were conducted us- the numbers of changes at different codon posi- 
mg the Branch and Bound option in PAUP for the tions. 
ì acıntosh computer (version 3.0L; Swofford, To evaluate further the relative robustness of 
990), which guarantees finding all most parsi- clades found in the most parsimonious trees, strict 
an trees. PAUP was also used in computing consensus trees were constructed of all trees up to 
ee trees (including “combinable compo- one step longer than the shortest trees, then of all 
ee trees: Bremer, 1990), optimizing character trees two steps longer, and so forth, until the con- 
i e changes (Swofford & Maddison, 1987), per- sensus tree eventually collapsed to an unresolved 
orming bootstrap analyses (100 replicates: Fel- bush. This procedure yields a “decay” index, in- 
eee 1985; Sanderson, 1989), differentially dicating the number of steps that must be added 
in ‘Siting characters (including use of a step matrix before each clade present in the minimum length 
cule Ee aenting the Tv: Ts weighting), and cal- trees is no longer unequivocally supported (cf. Bre- 
ie ing character and tree statistics, such as con- mer, 1988; Hillis & Dixon, 1989). 
ency index (CI; Kluge & Farris, 1969) and In order to explore more fully the structure in 
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FIGURE 1. Four most parsimonious trees of 202 steps each, rooted by Heuchera. The number associated with 
each internal branch indicates the number of unambiguous nucleotide changes along the branch. In 1A the pecs 
of transversions is shown in the first set of parentheses, and changes in first and second positions within codons ar 
indicated in the second set of parentheses: (1) = one first position change; (2) = one second position change. 


the data set, two other procedures were utilized. 
First, the length of the most parsimonious trees for 
the rbcL data was compared to the distribution of 
most parsimonious trees based on 50 “random” 
data sets generated under the null model of Archie 
(1989b). This was accomplished using the “shue” 
option in MacClade, which maintains the same 
number of states for each character observed in 
the real data set, but randomly assigns these states 
to the taxa. To save time, analysis of each of the 
50 random data sets was carried out using a heu- 
ristic search strategy in PAUP (‘‘closest” add se- 
quence, one tree held with each addition, TBR 
branch-swapping, and MULPARS). Second, PAUP 
was used to generate a histogram of tree lengths 
of 10,000 randomly sampled trees (there being too 
many taxa in the analysis to examine all trees). 
This allowed us to determine whether the observed 
distribution showed significantly greater left-hand 
skewness than tree-length distributions obtained 
from randomly generated data matrices (Hillis, 
1991; Huelsenbeck, 1991). The comparison is 
based on the g, skewness statistic (i.e., the third 


central moment divided by the cube of the standar d 
deviation; Sokal & Rohlf, 1981), for which critical 
values (P = 0.05) were provided by David Hilts 
(pers. comm.). RETE E 
Rooted trees were obtained by positioning ei 
root along the branch connecting Heuchera to : 
rest of the network (simultaneous resolution; Ma 3 
dison et al., 1984). Initially, we planned to position 
the root along the branch connecting the ua 
rosid taxa to the rest of the network, under í 
assumption that asterid taxa would join the mon 
along a single branch and that the root should 5 
be placed within Asteridae. However, it emerg h 
that the rosids were not attached by a single br pe 
to the asterid taxa in the analysis (see below). 
decision to root the tree along the Heuchera brane 
assumes only that Saxifragaceae is a more 
rosid clade than Cornaceae and/or Araliaceae. 


RESULTS 


$ ars l 
Parsimony analysis of the matrix in To in 
resulted in the four trees of 202 steps shown 
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Figure 1. The strict consensus of these trees (iden- 
tical to the combinable component tree) is shown 
in Figure 2. The CI of 0.515 (RI = 0.470) indicates 
that there are about twice as many character state 
changes as there would be in the absence of ho- 
moplasy. Although this is a rather low CI for 12 
taxa, there is still not a statistically significant de- 
parture from the expected value of 0.64 based on 
an analysis of 60 published data sets (Sanderson 
& Donoghue, 1989). 

Note that in Figure 1 we have indicated the 
number of character state changes that can be 
unambiguously assigned to each internal branch, 
and for one tree (Fig. 1A) we also show the number 
of unambiguous changes that are anything other 
than third-position transitions. As expected from 
the matrix, the vast majority of changes in the 
Most parsimonious trees occur at third positions. 
Separate analyses of characters at diffeent codon 
Positions confirm that the basic structure seen in 
the most parsimonious trees is primarily determined 
by these third position characters. When third po- 
sitions are analyzed by themselves, one tree of 163 
oh (CI = 0.546) is obtained, which is identical 
to Figure 1B. In contrast, when only the 14 first 
and second position characters are analyzed the 
consensus of the 55 trees of 31 steps (CI = 0.484) 
completely unresolved. Weighting first and sec- 
ond position changes twice as much as third position 
changes allows discrimination among the four trees 
chara l (A: 240 steps; B: 241; C: 238; D: 
), but we see no justification for such weighting 
In this case, 

i Also as expected, the number of transitions far 

xceeds the number of transversions on the most 
eee trees, yet the percentages of the two 
“th sof changes on the trees differ somewhat from 
Es expectation based on the matrix. From the 
a 16. it can be seen that there are a minimum 
as Be ssi and 28 transversions, which means 
et e might expect about 34 of the changes on 
rai ee to be transitions, However, the tree shown 
ieee 1A necessitates 102 unambiguous tran- 
e> be 45 unambiguous transversions. Thus, 
saciid ut 34 of the unambiguous changes are 
ease meaning that there are relatively fewer 
ph ‘ons, and relatively more transversions, than 

© might have expected based on the matrix. 
feck the different numbers of changes of the 
Pe tin differential weighting of transversions 

-4S ‘fansitions had little effect on the outcome. 
tr arms transversions a weight of 1.3 (see above) 
ká te trees, corresponding to Figures 1A, 1C, 
Was The consensus of these trees is identical 
ees consensus based on equal weighting (Fig. 2). 

other words, this weighting scheme has the effect 


Donoghue et al. 337 


rbcL Sequences of Dipsacales 


* 


SOPPPEEOLES E 


8-10 


FIGURE 2. Strict consensus of the four most parsi- 
monious trees in Figure 1. Numbers associated with the 
internal branches indicate the range of unambiguous char- 
acter changes, and the percentage of times that the branch 
was recovered in 100 bootstrap samples. Taxa marked 
by an asterisk are traditionally placed in Caprifoliaceae; 
taxa marked by (R) are assigned to Rosidae. 


of rendering one of the most parsimonious trees 
slightly less parsimonious (see Albert et al., 1992, 
for a similar example in Clarkia). Even a weighting 
of 2:1 had no effect on the consensus result. 

Along each branch of the consensus tree we 
indicate the range of unambiguous character 
changes and the percentage of times the branch 
was recovered in 100 bootstrap replicates. Only 
one branch not seen in the strict consensus ap- 
peared in over 50% of the bootstrap runs: Vale- 
riana and Dipsacus were directly united in 52% 
of the trees. Two clades are seen in over 95% of 
the samples: Lonicera plus Symphoricarpos (99%) 
and Sambucus plus Adoxa (97%). The connection 
of Lonicera and Symphoricarpos with Valeriana 
and Dipsacus is not far behind, being found in 
92% of the replicates. 

The order of “‘decay” of the six components 
present in the most parsimonious trees is shown in 
Figure 3. Two clades are lost when the six trees 
just one step longer are considered, and the re- 
maining four clades are maintained with the ad- 
dition of 21 trees of length 204 and 57 trees of 
205. When the 182 trees less than or equal to 
206 steps are considered, the connection between 
Viburnum and the Sambucus—Adoxa clade is lost. 
With the addition of the 142 trees of 207 steps 
the consensus tree no longer maintains the con- 
nection between Sambucus and Adoxa, although 
this link is still found in 99% of the 324 trees. 
Only the connection between Lonicera and Sym- 
phoricarpos is still present in the consensus of all 
594 trees of 208 steps or less, and all structure 
is lost at 209 steps, only 7 steps from the most 
parsimonious trees. 


, 
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beh? buble buble 


4 trees 10 trees 
202 steps < 203 steps 


182 trees 
< 206 steps 


Lett babel buble 


324 trees 
< 207 steps 


594 trees 
< 208 steps 


1021 trees 
< 209 steps 


FIGURE 3. Trees showing the order of “decay” of clades. Each tree is the strict consensus of the trees that ot 
less than or equal to the number of steps shown. The basal branch in all trees is an artifact of rooting the netwo 


along the branch connecting Heuchera. 


In Figure 4 the tree lengths of 50 randomized 
data sets are compared to the length (202) of the 
most parsimonious trees obtained for the rbcL data. 
Clearly, the real data set falls far (several standard 
deviations) outside the normal distribution for data 
sets randomized using the method suggested by 
Archie (1989b). Whether this result will hold for 
other null models is uncertain, but it may suggest 
that there is phylogenetic structure in the rbcL 
data. 


The distribution of the lengths of 10,000 ran- 


15 


WUMBER OF ANALYSES 


E 3 à i § Ei] g E $ 

TREE LENGTH 
FIGURE 4. Histogram of tree lengths showing the 
position of the most parsimonious trees (202 steps) for 
the rbcL data in relation to parsimonious trees for 50 
randomized data sets (sensu Archie, 1989b); arrow in- 


dicates the most parsimonious trees found with the real 
data set. 


k 


domly generated trees is shown in Figure 5. Note 
that the shortest random tree recovered was of 
length 220—18 steps longer than the shortest 
trees for the rbcL data. Although this means that 
none of the 1,021 trees known (from the decay 
experiments) to require 209 steps or less were 
actually found in the sample of random trees, it r 
presumed that 10,000 trees is sufficient to revè 

the shape of the underlying distribution for all pe 
(D. Hillis, pers. comm.). The skewness of this dis- 
tribution, measured by g,, is —0.78, which far 


400 


300 


NUMBER OF TREES 


200 


eR S e se Ss 
TREE LENGTH 
Ficure 5. Histogram of tree lengths for pe ran 
domly generated trees, showing the degree of 
skewness (g, = —0.78). 
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DIPSACALES 


LR 


203 steps 


CAPRIFOLIACEAE S.L. 


tl fle 


G 210 steps 


CAPRIFOLIACEAE 


LG IGE fee 


B eps 


ASTERIDAE 


ti thle 


Ww 


ob : aba Lengths of alternative trees, with the following topological constraints imposed: (A) Dipsacales mono- 
y rs 4 c; (B) Caprifoliaceae (Sambucus, Viburnum, Lonicera, Symphoricarpos) monophyletic; (C) Caprifoliaceae plus 
a monophyletic; (D) Asteridae monophyletic. When more than one tree was obtained, the strict consensus is 


shown 


Senge the P = 0.05 critical value for data sets 
W 10 to 15 taxa with 50 to 100 characters (D. 
» 1991, pers. comm.). Again, this appears to 


none significant phylogenetic signal in the rbcL 


Discussion 


DIPSACALES AND CAPRIFOLIACEAE SENS. LAT. 


a ae Dipsacales have generally been considered 
rt. aa ’ natural element within Asteri- 
Re egenitz, 1977, 1992). However, none of 
h Hea parsimonious trees support the mono- 
Phyly of Dipsacales (Fig. 1), even though the con- 
aesa of these trees (Fig. 2) seems to allow this 
Possibility. Instead, the clade consisting of Meny- 
Ey etoeto Asterales sens. lat.; Olmstead 
Rosid, 992) and Hedera (representing Apiales, 
Ree ae) is intercalated between two clades of dip- 
hie ean taxa, namely a Viburnum—Sambucus— 
h a clade and a line comprising Lonicera, Sym- 
Phoricarpos, Valeriana, and Dipsacus. 
regan at face value, this result supports the 
ey sed on immunology (Hillebrand & Fair- 
i sii 1970a, b) and morphology, that these 
ogh ines of Dipsacales are not closely related (Don- 
ue, 1983). It is also generally consistent with 


restriction site studies of the inverted repeat (Down- 
ie & Palmer, 1992) and other regions of the chlo- 
roplast genome (Sytsma, Smith, & Donoghue, un- 
published). However, it is important to recognize 
that the non-monophyly of Dipsacales is not a 
strongly supported result of our analysis. In fact, 
trees in which Dipsacales are monophyletic (Fig. 
6A) require only a single additional step. Note that 
in this case the Menyanthes-Hedera clade ap- 
pears as the sister group of Dipsacales, an arrange- 
ment that is more consistent with the broader anal- 
ysis of rbcL sequences of Olmstead et al. (1992). 

None of our most parsimonious trees are con- 
sistent with the monophyly of Caprifoliaceae in the 
traditional sense (that is, excluding Adoxa, Val- 
erianaceae, and Dipsacaceae), and this conclusion 
is very strongly supported in the present analysis. 
In particular, Viburnum and Sambucus are united 
with Adoxa, and this clade is separated from Lo- 
nicera and Symphoricarpos (of Caprifoliaceae sens. 
str.), which are united instead with Valeriana and 
Dipsacus. Trees in which Sambucus, Viburnum, 
Lonicera, and Symphoricarpos are united in a 
clade require a minimum of 216 steps, 14 more 
steps than the shortest trees (Fig. 6B). Even if 
Adoxa is considered part of Caprifoliaceae (as it 
sometimes has been; see Donoghue, 1983), it would 
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still be better to consider the latter to be paraphy- 
letic, because trees with this expanded Caprifoli- 
aceae forming a clade require a minimum of 210 
steps (Fig. 6C). 

Regarding the status of Caprifoliaceae our re- 
sults are congruent with morphological characters. 
Donoghue (1983) summarized earlier evidence, and 
research available since that time—for example, 
on pollen (Bohnke-Giitlein & Weberling, 1981; 
Donoghue, 1985), nectaries (Wagenitz & Laing, 
1984), and the tapetum (Weberling & Hilden- 
brand, 1986)—has strengthened the idea that there 
are independent and unrelated lines of Caprifoli- 
aceae. This view is also consistent with a phenetic 
analysis of morphological data (Hsu, 1983), with 
chromosome structure (Sax & Kribs, 1930), with 
secondary chemistry (Bohm & Glennie, 1971), and 
with immunological evidence (Hillebrand & Fair- 
brothers, 1970a, b). Finally, our results are con- 
gruent in this respect with the results of other 
molecular studies (Olmstead et al., 1992; Downie 
& Palmer, 1992). 


CAPRIFOLIACEAE SENS. STR., 
VALERIANACEAE, AND DIPSACACEAE 


The result that genera of Caprifoliaceae (in the 
traditional sense) do not form a monophyletic group 
is partly a function of the solid connection between 
Caprifoliaceae sens. str., Valerianaceae, and Dip- 
sacaceae. Indeed, the Lonicera—Symphoricar- 
pos~Valeriana—Dipsacus clade is one of the best 
supported in the entire analysis, judging by the 
number of unambiguous character changes, the 
bootstrap value of 92%, and the high “decay” 
index (collapsing in trees 6 steps longer than the 
shortest trees). This connection is also supported 
by many morphological features, including reduc- 
tion of carpellary vascular traces (Wilkinson, 1949), 
aspects of pollen morphology (Donoghue, 1985), 
and nectary structure (Wagenitz & Laing, 1984). 

Within this clade there is very strong support 
for a direct connection between Lonicera and Sym- 
phoricarpos. These genera are united by 8-10 
unambiguous changes in the most parsimonious 
trees, and they appeared together in 99% of the 
bootstrap samples. Furthermore, this link lasted 
the longest in the decay analyses, appearing in all 
594 trees less than or equal to 208 steps. At 209 
steps it is possible to find trees in which Sym- 
phoricarpos and Lonicera are dissociated, with 
either one united with Valerianaceae and Dipsa- 
caceae, 

The strong connection between Lonicera and 
Symphoricarpos is somewhat surprising. These 
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two genera are not unequivocally connected by 
morphological characters and they have generally 
been placed in separate tribes (see Wilkinson, 1949; 
Donoghue, 1983). Lonicera is usually allied with 
Leycestaria in the tribe Lonicereae. More impor- 
tantly, Symphoricarpos has often been placed in 
the tribe Linnaeeae (e.g., Fukuoka, 1972; Hara, 
1983), whose other members appear to be more 
closely related to Valerianaceae and Dipsacaceae 
than they are to other genera of Caprifoliaceae 
sens. str. (Donoghue, 1983, in prep.). If Sym- 
phoricarpos really belonged with Linnaeeae, and 
if Linnaeeae are more closely related to Valeriana- 
ceae and Dipsacaceae, then it should not have 
attached directly to Lonicera in our analysis, but 
instead to Valeriana and Dipsacus. Other molec- 
ular data seem to’ confirm that Symphoricarpos 
is more closely related to Lonicera than it is to 
genera of Linnaeeae. The two genera are united 
in the trees of Downie & Palmer (1992), and it 
appears that they may also be united (along with 
Triosteum) by the loss of the chloroplast gene clpP 
(Downie, Donoghue, Sytsma, & Palmer, unpub- 
lished). 

Our results for Lonicera and Symphoricarpos 
might be considered evidence of the monophyly of 
Caprifoliaceae sens. str., but this conclusion would 
be premature in view of the limited sample of taxa 
included in the analysis. As we noted above, one 
of the strongest results of morphological analyses 
is that Caprifoliaceae sens. str. are paraphyletic, 
with the tribe Linnaeeae more closely related to 
Valerianaceae and Dipsacaceae than to other Cap- 
rifoliaceae sens. str. (Donoghue, 1983, in prep.): 
Some of the strongest evidence is provided by flor 
anatomy, especially the abortion of carpels during 
the development of the ovary (Wilkinson, 1949). 
The link between Linnaeeae and Valerianaceae and 
Dipsacaceae is also evident in the restriction site 
studies of Downie & Palmer (1992), although m 
this case the one representative of ere 
(Kolkwitzia) is actually connected directly W 
Valeriana. It seems likely that if representatives 
of Linnaeeae were included they would be link 
with Valeriana and Dipsacus rather than W í 
Lonicera and Symphoricarpos. Sequencing y 
representatives of Linnaeeae (Dipelta, Kolkwit2® 
Linnaea) is in progress. A 

In our analysis the position of Valeriana is ™™- 
certain. In three of the four most parsimon! 
trees (Fig. 1) it is united directly with Dipsacus 
but in the fourth it is linked with Lonicera a 
Symphoricarpos. Morphological data support s 
close connection between Valerianaceae 


sacaceae (e.g., Hofmann & Göttmann, 1990; Don- 
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oghue, in prep.), as does the more inclusive rbcL 
analysis of Olmstead et al. (1992), which again 
suggests that there may be a problem of limited 
sampling of taxa. Inclusion of additional Caprifo- 
liaceae sens. str. might strengthen the connection 
between Valerianaceae and Dipsacaceae. 


VIBURNUM, SAMBUCUS, AND ADOXA 


Viburnum is linked with Sambucus plus Adoxa 
in all of our most parsimonious trees, but this 
connection is not especially srong: 2—4 unambig- 
uous character changes mark this branch, it ap- 
peared in 77% of the bootstrap trees, and it de- 
cayed at 206 steps. Morphological data also tend 
to unite Viburnum with Sambucus and Adoxa 
(under most outgroup arrangements), but by itself 
the evidence does not inspire great confidence 
(Donoghue, 1983, in prep.). It is certainly en- 
couraging that the molecular and morphological 
data point to the same conclusion. 
The connection between Adoxa and Sambucus 
'S quite robust. There are 6-7 unambiguous char- 
oa changes below this node, and it appeared in 
1% of the bootstrap samples. In the decay anal- 
ais this connection lasted until trees of 207 steps 
ne added, making it the third most robust result 
re toa Curiously, this clade decayed before 
at teri a—Symphoricarpos—Valeriana—Dip- 
tee ee e, even though the latter had a somewhat 
tstrap value. 
oe to Caprifoliaceae (e.g., Eichler, 1875; 
iis ah 1873), Adoxa has been allied with Sax- 
a Ps (Hutchinson, 1973; Warming & Mo- 
val 9), Araliaceae (de Candolle, 1830), and 
wks r (Dahlgren, 1975), although it has gen- 
Ware een placed in its own family (reviewed in 
be ee nitz, 1977; Donoghue, 1983). The rbcL data 
des against these other connections. Thus, the 
wii = A in which Adoxa is linked directly 
iban ik era require 213 steps (11 steps longer 
Pes most parsimonious trees), and it takes 214 
x connect Adoxa with Heuchera. 
the, authors have noted similarities between 
hee om Sambucus (e.g., Eichler, 1875; Hal- 
RP ; Takhtajan, 1980), and this connection 
sai Da supported by preliminary cladistic anal- 
ical E nA AR 1983, in prep.). One morpholog- 
ie ae in particular has seemed compelling, 
dena the presence of Adoxa-type embryo sac 
dition Pment in Adoxa and Sambucus. This con- 
ny extremely rare in angiosperms and is oth- 
haa very widely scattered in occurrence (Ma- 
waives 1946; Donoghue, 1983). Its distribution 
ambucus is still somewhat unclear, and it 


is possible that it will be discovered in some Vi- 
burnum species, inasmuch as several species are 
reported to have bisporic development (Sunesson, 
1933). Other possible indicators of a close rela- 
tionship include compound leaves (depending on 
outgroups) and extrorse anthers (although there is 
some variation in Sambucus and the condition in 
Adoxa is difficult to score). Nevertheless, most 
authors (e.g., Sprague, 1927) have hesitated to 
connect the two genera, largely because Adoxa is 
so different (presumably highly derived) in habit 
and in floral morphology. However, Thorne (1983) 
adopted Donoghue’s (1983) suggestion that Adoxa 
and Sambucus (and possibly Viburnum) be treated 
as the family Adoxaceae. 

One outstanding problem concerns the relation- 
ships of recent additions to Adoxaceae (Wu, 1981; 
Wu et al., 1981; Hara, 1981, 1983). Although 
Tetradoxa omeiensis (=~Adoxa omeiensis) is very 
similar to Adoxa moschatellina, and is probably 
closely related, Sinadoxa differs in many respects 
and may be more closely related to Araliales than 


to Adoxa. 


ASTERIDAE AND ROSIDAE 


Conclusions regarding the broader relationships 
of dipsacalean plants are tenuous, because much 
of the relevant structure seen in the most parsi- 
monious trees disappears when slightly longer trees 
are taken into consideration. Nevertheless, it is 
interesting to note that the most parsimonious trees 
support a clade including the dipsacalean taxa, 
along with Hedera plus Menyanthes. The former 
is a placeholder for Araliales, and the latter rep- 
resents an Asterales sens. lat. clade (comprising 
Asteraceae, Menyanthaceae, Goodeniaceae, Cal- 
yceraceae, Campanulaceae, and Lobeliaceae), which 
is consistently supported in the more inclusive rbcL 
analyses of Olmstead et al. (1992). Although the 
exact arrangement of these groups in relation to 
dipsacalean taxa is uncertain, the present analysis 
is consistent with the view that Asterales sens. lat. 
and Araliales are closely related, and that together 
they are related either to the Viburnum—Sambu- 
cus-Adoxa clade or to the Caprifoliaceae sens. 
str.—-Valerianaceae—Dipsacaceae clade. 

This result is especially significant in light of the 
inclusion of Spigelia (Loganiaceae) in the analysis, 
which is intended to represent the remainder of 
the Asteridae, based on the broader analyses of 
Olmstead et al. (1992). If Asteridae were mono- 
phyletic, Spigelia would be linked with Meny- 
anthes and the dipsacalean taxa, to the exclusion 
of the rosids. The shortest trees of this type require 
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206 steps, four steps more than our most parsi- 
monious trees (Fig. 6D). Instead, Cornus and Hed- 
era are variously nested among Asteridae in our 
most parsimonious trees (Fig. 1), which implies that 
the Asteridae, as they are generally constituted, 
are not monophyletic. This conclusion is consistent 
with the larger analyses of rbcL (Olmstead et al., 
1992). 

Although the standard view of morphological 
systematists has been that Asteridae form a natural 
group, Takhtajan (1987) proposed that it be sub- 
divided into two major lines: Asteridae sens. str. 
(Asterales sens. lat. and Dipsacales) and Lamiidae 
(the remaining Asteridae). These lines, he believes, 
originated independently from separate lines of 
Rosidae. Phylogenetic analyses of the morpholog- 
ical evidence are needed before molecular and mor- 
phological results can be compared critically. 

In the meantime, it is interesting to note that 
while the results of the present study are consistent 
with many aspects of other molecular analyses, 
they also differ in several ways. In trees presented 
by Olmstead et al. (1992) Araliales are linked di- 
rectly with Dipsacales, rather than with Asterales 
sens. lat. Similarly, the trees of Downie & Palmer 
(1992) show Araliales connected to Caprifoliaceae 
sens. str., Valerianaceae, and Dipsacaceae. Al- 
though these studies include more taxa outside of 
the groups of immediate concern here, they also 
contain fewer potentially close relatives. The res- 
olution of this issue requires a more complete sam- 
ple of relevant taxa as well as additional molecular 
and morphological evidence. 

It is important to recognize that all three of the 
molecular studies suggest that Asteridae in their 
traditional sense are not monophyletic, and they 
imply that some “‘asterid” features such as sym- 
petaly may have evolved independently in Asterales 
sens. lat., Dipsacales, and Lamiidae, or may have 
been lost in Araliales and possibly Cornales. Here, 
the ontogenetic studies of Erbar (1991) are of 
special interest. She has demonstrated a difference 
between “‘early sympetaly” (corolla tube initiated 
before petal primordia) and “‘late sympetaly”’ (co- 
rolla lobes initiated as separate primordia and con- 
nected later), which largely corresponds with major 
clades of Asteridae implied by molecular data. In 
particular, early sympetaly is found in Asterales 
sens. lat. (including Campanulales and Menyan- 
thaceae) and Dipsacales, while late sympetaly oc- 
curs in virtually all other Asteridae (i.e., in La- 
miidae sensu Takhtajan, 1987). Rubiaceae, 
Oleaceae, and several Acanthaceae are exceptions 
to this correspondence in having early sympetaly. 
Unfortunately, the development of sympetaly has 
not yet been studied in Viburnum, Sambucus, or 


Adoxa. These observations are consistent with (but 
do not necessitate) the idea that sympetaly evolved 
independently at least twice. 

Erbar (1991) also noted evidence (see Erbar & 
Leins, 1988) that in some Araliales the corolla is 
initiated as a ring primordium. Because this ring 
does not develop into a conspicuous tube, however, 
the mature flowers of these plants appear to have 
separate petals. This similarity to early sympetaly 
can be interpreted in several ways. It may be that 
the condition in Araliales is a precursor of early 
sympetaly in Asterales sens. lat. and Dipsacales. 
Alternatively, the apparently separate petals of 
Araliales may have been derived from the early 
sympetalous condition (i.e., sympetaly was “lost”. 
The latter hypothesis will be favored if phylogenetic 
studies continue to nest Araliales among sympet- 
alous groups (see Olmstead et al., 1992). 


CONCLUSIONS 


The analyses reported here suggest a number 
of changes in the circumscription of taxa. First, it 
is clear that the genera of Caprifoliaceae do not 
form a monophyletic group, and the traditional 
concept of this family must be abandoned if we 
hope to achieve a truly phylogenetic system of 
angiosperms (Donoghue & Cantino, 1988). In par- 
ticular, the genera of Caprifoliaceae sens. str. are 
more closely related to Valerianaceae and Dipsa- 
caceae than they are to Viburnum and Sambucus, 
and Sambucus (very probably) and Viburnum (with 
less certainty) are more closely related to Adoxa. 
The smallest monophyletic group that would in- 
clude the genera traditionally assigned to Capri- 
foliaceae would also include the rest of the Dip- 
sacales (Valerianaceae, Dipsacaceae, and 
Adoxaceae), and possibly also the Asterales sens. 
lat. and Araliales. 

The connection between Caprifoliaceae sens. St 
Valerianaceae, and Dipsacaceae should be formally 
represented, but a decision on the best way ae 
this must await the publication of the morpholog'¢ 
results and additional molecular analyses. The r be 
sequence data are still too limited to evaluate pr op: 
erly the morphological hypothesis that Caprifolia- 
ceae sens. str. are paraphyletic, with the t 
Linnaeeae more closely related to Valerianacea® 
and Dipsacaceae. If this is borne out by further 
studies, it may be best to merge Caprifoliacea® 
sens. str., Valerianaceae, and Dipsacaceae into & 
single family. 

Although the status of Caprifoliaceae sens. we 
is still uncertain, Lonicera and Symphoricarpe 
seem to be more closely related than one 4 
have expected based on previous taxonomic trea 
ments of the family (see discussion above). It sh 
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be noted, however, that because ovary develop- 
ment in Symphoricarpos does not correspond ex- 
actly with Lonicereae or with Linnaeeae (Wilkin- 
son, 1949), its exact relationships based on 
morphology are still uncertain (Donoghue, in prep.). 
Our analyses strongly suggest that rbcL sequences 
of representatives of additional genera of Capri- 
foliaceae sens. str. would be very valuable and 
might result in tribal realignments. 

The close relationship between Sambucus and 
Adoxa should be represented in the phylogenetic 
system, perhaps by expanding Adoxaceae. This has 
been suggested based on morphological evidence 
(Donoghue, 1983, in prep.), but these data have 
not been entirely convincing, owing especially to 
the extreme modification/reduction of Adoxa. The 
rbcL data provide strong support for this connec- 
tion and argue against alternative hypotheses that 
Adoxa may be related to Araliales or Saxifraga- 
ceae, Recently described Adoxaceae, especially 
Sinadoxa, must be evaluated, as should additional 
representatives of Sambucus. This phylogenetic 
information will be especially critical in understand- 
ing the nature of the genetic and developmental 
changes associated with morphological shifts in the 
evolution of Adoxa-like plants. 

The position of Viburnum is still somewhat un- 
certain. It is connected with Sambucus and Adoxa 
in our most parsimonious trees, and there is also 
some morphological support for this arrangement. 
Although the evidence is not overwhelming in either 
case, the congruence is satisfying. Viburnum might 

placed in Adoxaceae along with Sambucus and 
Adoxa, but until its relationships are more confi- 
dently resolved there may be some merit in placing 
Viburnum in its own family, the Viburnaceae (see 
Dahlgren, 1980). In the end, the assignment of 
rank is arbitrary. 

_ Our results are encouraging as regards the qual- 
my: of rbcL data. Although the consistency index 
'S quite low for a data set of this size, randomized 
data sets and tree-length distributions both suggest 
that considerable phylogenetic information is pres- 
ent in the data. The bootstrap and decay analyses 
Suggest that this signal may reflect the presence 
of several very well supported clades: especially 
the connections between Lonicera and Symphori- 
“arpos, between these genera and Valeriana and 

'psacus, and between Sambucus and Adoxa. 
Most of the rest of the connections are not strongly 
Supported and decay with the addition of only a 
Single step. Overall, these results are especially 
encouraging with respect to potential resolution of 
the relationships of dipsacalean taxa using rbcL 
Sequences. 

Our analyses also bear on the general issue of 


differentially weighting characters or character state 
changes. Our main phylogenetic results are sup- 
ported by changes at third positions within codons; 
there are few changes at first and second positions 
and these show considerable conflict with one an- 
other. We see no reason to disregard or even to 
downweight third sites in this case, and we caution 
against a priori judgements of this sort. Likewise, 
we found that much of the phylogenetic structure 
resided in the much more numerous transition sub- 
stitutions, and saw no evidence that transitions were 
less reliable than transversions. 

Comparison of this study with the more inclusive 
rbcL analysis of Olmstead et al. (1992) raises im- 
portant issues concerning the effects of the sam- 
pling of taxa in phylogenetic studies. The smaller 
sample considered here allows a much more de- 
tailed examination of the structure in the data set. 
Exact parsimony solutions can be obtained using 
branch-and-bound algorithms, which also ensure 
an accurate picture of the rate of decay of clades, 
and so on. Studies using randomized data sets and 
tree-length distributions are also much easier to 
carry out with smaller data sets. Some of these 
sorts of analyses can be extended to larger analyses, 
and we hope that exploration of the robustness of 
clades (beyond the use of the bootstrap) will soon 
become routine. In the meantime, it is difficult to 
ascertain the reliability of the results of the larger 
analyses, and one must be especially cautious in 
drawing evolutionary conclusions based on such 
studies (see Maddison, 1991). 

On the other hand, smaller studies are subject 
to error by virtue of leaving out relevant taxa 
(Donoghue et al., 1989; Olmstead et al., 1992). 
In particular, in any analysis of a small portion of 
Asteridae (for example) one may not discover glob- 
ally more parsimonious solutions, which might be- 
come apparent with the addition of taxa outside 
the group of interest. Of course, this warning ap- 
plies to larger studies as well, as does the problem 
of leaving out taxa within the group of interest that 
might have an important impact on the outcome. 
If the sample of taxa within the group is too sparse, 
errors might arise by virtue of highly divergent 
branches “attracting” one another due to the chance 
accumulation of mutations (see Felsenstein, 1978). 
Increasing the sample of taxa might effectively 
“shorten” such long branches, thereby decreasing 
the likelihood of such errors. This sort of consid- 
eration could perhaps be translated into a scheme 
for selecting a set of taxa in the first place. 

In our analysis there is a better representation 
of critical dipsacalean taxa, which probably allows 
our study to reflect more accurately relationships 
among Dipsacales than the broader analysis of Olm- 
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stead et al. (1992). On the other hand, our rep- 
resentation of other major lines is much poorer 
than in the analysis of Olmstead et al. (1992). Our 
failure to achieve any clearer picture of the broader 
relationships may be largely a function of the lim- 
ited sample of taxa aside from Dipsacales, rather 
than any fundamental problem with the use of rbcL 
sequences in addressing the more inclusive phy- 
logenetic questions. 

The congruence between molecular and mor- 
phological evidence on the phylogeny of Dipsacales 
is gratifying, if not very surprising (Donoghue & 
Sanderson, 1992). As is often the case, it is clear 
that the different data sets provide consistent re- 
sults, but that the results show different levels of 
resolution. For example, in our rbcL analysis the 
position of Valeriana is unresolved, whereas sev- 
eral morphological characters link Valerianaceae 
with Dipsacaceae. On the other hand, rbcL se- 
quences strongly connect Symphoricarpos with 
Lonicera, whereas the relationships of these genera 
are not clearly resolved in morphological analyses 
(Donoghue, in prep.). Under these circumstances 
it will be especialy valuable to combine data sets, 
so that all of the evidence can be considered at 
once (Barrett et al., 1991; Donoghue & Sanderson, 
1992). 
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